extern int trapru_(E_fp f, doublereal *a, doublereal *b, integer *m, doublereal *trule);
extern int xtrapru_(E_fp f, doublereal *a, doublereal *b, integer *m, doublereal *trule);
